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ABSTRACT 

A typical semiconductor-based optoelectronic device, such as a diode laser, consists of 
three subsystems: an optical field, an electron-hole plasma (EHP), and a host crystal lat- 
tice. The physics of such a device involves the interplay of optical, electrical and thermal 
processes. A proper description of such a device requires that all three proeesses are treated 
on equal footing and in a self-consistent fashion. Furthermore, since a semiconductor laser 
has intrinsic spatial inhomogeneity, such a self-consistency naturally leads to a set of partial 
differential equations in space and time. 

There is a significant lacking of research interest and results on the transport aspects of 
optical devices in the literature with only a few exceptions[l-3]. Even the most important 
carrier diffusion coefficient has not been properly derived and studied so far for optically 
excited plasma, while most of the work adopted results from electronics community where 
heavily doped semiconductors with mainly one type of carriers are dealt with. The corre- 
sponding transport equation for plasma energy or temperature has received even less atten- 
tion. In this talk we describe our recent results [4-6] on such a self-consistent derivation of 
temperature and carrier-density diffusion equations coupled with the lasing process. 

Starting from the microscopic semiconductor Bloch equations (SBEs) including the Boltz- 
mann transport terms in the distribution function equations for electrons and holes, we 
derived a closed set of diffusion equations for carrier densities and temperatures with self- 
consistent coupling to Maxwell’s equation and to an effective optical polarization equation. 
The coherent many-body effects are included within the screened Hartree-Fock approxima- 
tion, while scatterings are treated within the second Born approximation including both 
the in- and out-scatterings. Microscopic expressions for electron-hole (e-h) and carner-LO 
(c-LO) phonon scatterings are directly used to derive the momentum and energy relaxation 
rates. These rates expressed as functions of temperatures and densities lead to microscopic 
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expressions for self- and mutual-diffusion coefficients in the coupled density-temperature dif- 
fusion equations. Approximations for reducing the general two-component description of the 
electron-hole plasma (EHP) to a single-component one are discussed. In particular, we show 
that a special single-component reduction is possible when e-h scattering dominates over 
c-LO phonon scattering. The ambipolar diffusion approximation is also discussed and we 
show that the ambipolar diffusion coefficients are independent of e-h scattering, even though 
the diffusion coefficients of individual components depend sensitively on the e-h scattering 
rates. Our discussions lead to new perspectives into the roles played in the single-component 
reduction by the electron-hole correlation in momentum space induced by scatterings and 
the electron-hole correlation in real space via internal static electrical field. Finally, the 
theory is completed by coupling the diffusion equations to the lattice temperature equation 
and to the effective optical polarization which in turn couples to the laser field. 

The equations derived above are implemented in various limiting cases to a typical diode 
laser to study the consequences of nonlinear diffusion and the cross diffusion terms on laser 
behavior, especially the dynamic behavior of a diode laser under modulation. Detailed 
results will be presented by comparing with the standard rate equation results. 
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OUTLINE 

^ Introduction 

II. Equations of Motion 

III. Moment Equations 


n. Equations of Motion: 

Density Matrix Equations for Inhomog. Semiconductor (Kuhn, 1994 ) 


• Hamiltonian for a 2- band 2D e-h piasma 
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• The Heisenberg Equation for the bilinear op erato r combinations 
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Density Matrix Equations (continued) 


Non-diagonal density 


matrix elements in momentum representation 
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Non-diagonal density matrix elements in 'mixed' (Wigner) 
representation [integro -differentia I equations, (Kuhn, 1994)] 

K=(k+k')/2; q = k- k' 

Fourier transform w.r.t the relative momentum 
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Ames Research Center numerics 
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